Fast and robust two— qubit gates for scalable ion trap quantum computing 
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We propose a new concept for a two-qubit gate operating on a pair of trapped ions based on 
laser coherent control techniques. The gate is insensitive to the temperature of the ions, works also 
outside the Lamb-Dicke regime, requires no individual addressing by lasers, and can be orders of 
magnitude faster than the trap period. 
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Trapped ions constitute one of the most promising sys- 
tems to implement scalable quantum computation. 
In an ion trap quantum computer qubits are stored in 
long-lived internal atomic states. A universal set of sin- 
gle and two qubit gates is obtained by manipulating 
the internal states with lasers, and entangling the ions 
via the motional states 0- During the last years a re- 
markable experimental progress in building an ion trap 
uantum computer has allowed to realize two-qubit gates 
100 and also to prepare entangled states 0, 0, 0] . 
The ultimate challenge is now the development of scal- 
able ion trap quantum computing. Scalability is based 
on storing a set of ions, and moving ions independently, 
in particular to b ring to gether pairs of ions to perform 
a two-qubit gate [lfj, . Basic steps towards this goal 
have already been demonstrated experimentally [l^ . 

An important question to be addressed is to iden- 
tify the current limitations of the two-qubit gates with 
trapped ions (given the fact that one-qubit gates are sig- 
nificantly simpler with those systems) . The ideal scheme 
should: (i) be independent of temperature (so that one 
does not need to cool the ions to their ground state af- 
ter they are moved to or from their storage area); (ii) 
require no addressability (to allow the ions to be as close 
as possible during the gate so as to strengthen their inter- 
action) , and (iii) be fast (in order to minimize the effects 
of decoherence during the gate, and to speed up the com- 
putation). This last property has been identified [lj as 
a key limitation: in essentially all schemes suggested so 
far 0, El 0, 0, 0, E| one has to resolve spectroscopi- 
cally the motional sidebands of the ions with the exciting 
laser, which limits the laser intensity and therefore the 
gate time. 

The two-qubit gate between pairs of ions analyzed be- 
low solves the problem of speed by using mechanical ef- 
fects instead of spectral methods to couple the motion 
and internal states of the ions. In this way the new lim- 
its on the time of the quantum gate are those of laser 
control, which can be orders of magnitudes faster than 
the present limits dictated by trap design. Thus, the 
method proposed here is a significant step forward to- 
wards fast and efficient scalable quantum computations 
with trapped ions. 

Below we will first study the dynamics of a pair of 




FIG. 1: a) Trajectory in phase space of the center-of-mass 
state of the ion (X c ,Pc) (where (X c + iP c )/V2 = (a)) dur- 
ing the 2-qubit gate (solid line), connecting the initial state 
(black filled circle) to the final state (grey filled circle) at 
the gate time T. The time evolution consists of a sequence 
of kicks (vertical displacements), which are interspersed with 
free harmonic oscillator evolution (motion along the arcs). A 
pulse sequence satisfying the commensurability condition @ 
guarantees that the final phase space point is restored to the 
one corresponding to a free harmonic evolution (dashed cir- 
cle). The particular pulse sequence plotted corresponds to a 
four pulse sequence given in the text (Protocol I). Figure b) 
shows how the laser pulses (bars) distribute in time for this 
scheme. 



trapped ions under the influence of short laser pulses 
with varying directions. We will prove that there exist 
certain laser pulse sequences which perform a phase gate 
on the two qubits, while leaving the motional state un- 
changed. We illustrate this with two protocols for laser 
pulses: (i) a sequence of four pulses which gives a gate 
time of T — 1.08/v with v the trap frequency, and (ii) 
a protocol which allows us to perform a gate in a time 



T 



N p 2 ^ 3 jv where Np is the number of laser pulses. Fi- 
nally, we will complement our study of the gate dynamics 
with an analysis of possible errors, which includes fluctu- 
ations of the intensity or the duration of the pulses, and 
temperature. The gate will be shown to be extremely 
robust to these perturbations. 

We consider two ions in a one-dimensional harmonic 
trap, interacting with a laser beam on resonance. The 
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Hamiltonian describing this situational can be written 
as H = Hq + Hi, where Hq = u c a' a + v r b' b describes the 
motion in the trap and 

^]_ a + e tVc(a^ +a) + ^ Vr (b i ' +b) + 



Hi = 



^+ JVc(a 1 +a)-hv r (b f +b) 

— — — (7 o o z 

2 2 
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Here, v c = v and z^ r = \T&v c are the frequencies of the 
center of mass and stretching mode, respectively; a and 
b are the corresponding annihilation operators, and r\ c = 
77/ -s/2 and r\ r = r\\J 4/3 are proportional to the Lamb- 
Dicke parameter, r\. Note that the Rabi frequency 
is the same for both ions, since we have not assumed 
individual addressing. Also notice that replacing rj with 
is equivalent to reversing the direction of the laser 
beam. 

In the following we will consider two different kind of 
processes: (i) Free evolution, where the laser is switched 
off (Q = 0) for a certain time; (ii) Sequences of pairs of 
very fast laser pulses, each of them coming from opposite 
sides. If we denote by St the duration of a pulse and by 
Q the corresponding Rabi frequency, we are interested in 
the limit St — > with il St = it. Processes (i) and (ii) 
will be alternated (See Fig. [TJi): at time t\ a sequence 
of z\ pulses is applied, followed by free evolution until at 
time 7J2 another sequence of zi pulses is applied followed 
by free evolution and so on. The z k are integer numbers, 
whose sign indicates the direction of the laser pulses. 

For a pulse sequence, consisting of kicks interspersed 
with free harmonic time evolution (Fig. we write 
the evolution operator as U — U c U r , where U cr = 



n 



U Cjr (Atk, Zk) has contributions of the center-of- 



mass and relative motions, 



U c (tk,z k ) 

U r {tk,Z k ) 



-iz k ri r (b+b^ )((j\-oi i ) —iv r Atfcb+b 



The integers z k indicate the direction of the initial pulse 
in the sequence of pairs of very fast laser pulses, each of 
them coming from opposite sites. 

In order to fully characterize U, we only have to in- 
vestigate its action on states of the form |i)i|j)2|a) c |/3) r , 
where i,j = 0, 1 denote the computational basis, and \oi) 
and |/3) are coherent states. This task can be easily car- 
ried out once we know the action of U = Y\k=i U{<fik,Pk) 
on an arbitrary coherent state \a), where 



U(<f>k,Pk) 
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We obtain U\a) = e^la), where 
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FIG. 2: (a) Number of pairs of pulses , and (b) relative angle 
of the two laser beams required to produce a phase gate using 
the first exact scheme developed in the paper. 



with 9 k = J2 m =i <fim- 

The crucial point is to realize that if X^itLi Pk£ l6k = 
the motional state \a) after the evolution is the same as 
if there was only free evolution [Fig. ^] , and a global 
phase £ appears which does not depend on the motional 
state. Translating this result to the operators U c \a) and 
U r \(3), we obtain the conditions 
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If these commensurability conditions are satisfied, the 
motional state of the ion will not depend on the qubits 
and the evolution operator will be given by 

The value T is the total time required by the gate and 
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where Atkm = tk — t m . Therefore, if Eqs. © are ful- 
filled, and = 7r/4 we will have a controlled-phase gate 
(which is equivalent to a controlled-NOT gate up to local 
operations) which is completely independent of the initial 
motional state, i.e. there are no temperature require- 
ments. 

It is straightforward to show that for any value of 
the time T it is always possible to find a sequence of 
laser pulses which implements the gate, and therefore 
the gate operation can be, in principle, arbitrarily fast. 
The search for a sequence of pulses may be done numer- 
ically, or even semianalytically. In the following we give 
two simple protocols which are not optimized in order to 
reduce the number of pulses. 

The first protocol (Protocol I) requires the least num- 
ber of pulses to produce the gate in a fixed time T ~ 
1.08(27r/^). The recipe is illustrated in Fig. ^ which 
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FIG. 3: (a) Log-log plot of the number of pairs of pulses re- 
quired to produce a phase gate using protocol II, as a function 
of the duration of the gate, T, for a realistic value 4] of the 
Lamb-Dicke parameter, rj = 0.178. We plot both the exact 
result (solid line) and a rough estimate Np = 4Q(vT/2tt)~ 3 ^ 2 
(dashed line) based on perturbative calculations, (b) Max- 
imum relative displacement, X r (solid), and maximum mo- 
mentum acquired, P r (dashed line), for scheme II. These 
quantities are adimensionalized versions of the real observ- 
ables, X r = max[{x r (t)) /do], and P r max[(p r (f)} ao/h]. 



provides the phase space plots for the evolution of the 
motional state. The sequence of pulses is defined as 
(z n /N,t n ) = {( 7 ,-r 1 ),(l,-r 2 ),(-l,r 2 ),(-7,r 1 )}. Here 
< 7 = cos(#) < 1.0 is a real number, which may be 
introduced by tilting both lasers a small angle 6 with 
respect to the axis of the trap, so that no transverse mo- 
tion is excited. It is always possible to find a solution to 
Eq. © with n ~ 0.538(4) (2tt/z/) > t 2 > 0. The results 
for the performance of the gate are summarized in Fig. [21 
As shown in Fig. HJa), for realistic values of the Lamb- 
Dicke parameter |4j we only need to apply the sequence 
of pulses one or two times to implement a phase gate. 

The second protocol (Protocol II) performs the 
gate in an arbitrarily short time T. The pulses 
are now distributed according to (z n /N,t n ) = 
{(-2, -n), (3, -72,), (-2, -t 3 ), (2, 73), (-3, r 2 ), (2, n )}. 
The whole process takes a time T = 2t\ and requires 
N p = J2\ z n\ = 14-/V pairs of pulses. As Fig. [3] shows, 
the number of pulses increases with decreasing time as 
N p oc T~ 3 / 2 . 

In order to study the main potential limitations of our 
scheme, we define the error of the gate E in terms of the 

gate fidelity as E = l-Tr mot jQmotPmotQLrt} • Here 
Tr mot and Tr; n t denote traces over motional and internal 
degrees of freedom, and Q mot = Tr int |w(7r/4)J7 1 t cal | , 

depends on £/ rca i, the gate performed in the presence of 
imperfections. 

We now turn to a discussion of the possible sources 
of errors. A limiting factor for the gate is the anhar- 
monicities of the restoring forces. The more pulses we 
apply, the larger the relative displacement of the ions, 
as Fig. 02b) shows. When the ions become too close to 
each other, the increasing intensity of the Coulomb force 
can lead to a breakdown of the harmonic approximation 
which is implicit in Eq. l|T]l. In order to analyze this ef- 



fect, we have made a perturbation analysis for vT <C 1 
and found that such an anharmonicity causes an error 
E ~ |0.4ao/(i| 2 /(2m>T), where ao is the ground state 
size of the external potential and d is the ion separation 
in equilibrium. For typical parameters and imposing an 
error E ~ 10 -4 we obtain vT ~ 10~ 3 . A similar analy- 
sis could be applied to study anharmonicities of the trap 
itself. 

In addition, we have studied the influence of errors in 
the laser pulses of our scheme. Up to now, our analytical 
calculations assumed that the intensity of the laser is 
very large during each pulse, and that therefore one may 
neglect the influence of the trap during this process. To 
validate this assumption we have simulated numerically 
a system of two ions with only one vibrational mode. 
We have used the exact sequences developed above to 
produce the phase gate using only eight laser pulses. In 
Fig. 0Ja) we plot the error of the gate as a function of 
the duration of the laser pulse, r = n/2Cl. The longer 
the pulse, the more important the effect of the trap, and 
the larger the error. But even for relatively long pulses, 
we obtain a fidelity which is comparable to the results 
obtained in current setups 0, H, El • We have also studied 
the influence of noise in the intensity of the laser pulses, 
or, what is equivalent, random errors in its duration. The 
larger the amplitude of the error the lower the fidelity of 
the gate, as Fig. 0fb) shows. 

As mentioned before, the scheme is insensitive to tem- 
perature. If the commensurability condition Q is not 
perfectly satisfied due to, for example, errors in timing 
of laser pulses, or misalignment of the lasers, then the 
corresponding contribution to the gate error is 



with 



E = (Cf + C 4 + dC 2 - 6)/8 

Ci = cxp[-(l/2 + k b T/hv c )\2r] c C c \ 2 
C 2 = exp[-(l/2 + k b T/hu r )\ Vr C r \ 2 ] 



(6) 
(7) 



which is a smooth function of temperature T. 

Finally, we would like to make some remarks regarding 
the experimental implementation of this scheme. First, it 
is not necessary to kick the atoms using pairs of counter- 
propagating laser beams. The same effect (i.e. a change 
of sign in 77) may also be achieved in current experiments 
by reverting the internal state of both ions simultane- 
ously. One then only needs a laser beam aligned with 
the trap to kick the atoms, and another laser orthogonal 
to the axis of the trap to produce the NOT-gate. 

The second and more important remark is that it is 
possible to avoid errors in the laser pulses by using more 
sophisticated kicking methods. One possibility consists 
in using STIRAP [ljIHJ]. Only one of the qubit states 
would be connected by two on-resonance laser beams 
to a third atomic state, \e'). In the first part of the 
kicking process, the Rabi frequencies of the lasers fl a 
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FIG. 4: For the sequences of four kicks designed in the paper, 
usingQ rj = 0.178 and a duration of the gate T = 0(l/u), we 
have computed the dependence of the gate error with respect 
to the characteristics of the laser pulses. In (a) the error is 
plotted versus the duration of the pulse, r, rescaled using 
the period of the trap. It perfectly fits the estimate E — 
2(tu/2tv) 2 . In (b) the mean error of the gate is plotted against 
the amplitude of random errors in the duration of the laser 
pulses. Our errors are modeled setting tu = 7r/(2fi)(l + er/s), 
with random numbers uniformly distributed in [—1/2, 1/2]. 
In solid lines we show a visual aid which follows the formula 
E = 4e 2 . 



laser excitations tolerant against laser imperfections. In 
any case, the rapid experimental progress in laser control 
with very short pulses indicates that it may be possible 
soon to perform quantum gates with a very high speed. 
In addition, our scheme is independent of the tempera- 
ture, requires no addressability, and works beyond the 
Lamb-Dicke regime. 
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and fib are adiabatically switched on an off respectively. 
The momenta of both laser beams should be different, so 
that as we slowly proceed from Sl a /Qb — to the opposite 
regime f2b/0 Q ~ 0, the ions in the state |1) are completely 
transferred to the new dark state |e') and get a kick |1) — > 
e i(k a -kb)x | e ^_ ]\f e xt we must change the sense of the laser 
beams (k a ^ — > —k a ,b)i an d perform the adiabatic transfer 
from Q,b/Vl a ~ to £l a /^fc — 0. The total transformation 
will be |0) -> |0), |1) -> e 2i(fco-fci,)£|i). The advantages 
of this method are: (i) the system remains all the time 
in a dark state, avoiding spontaneous emission; (ii) the 
process is insensitive to fluctuations of the intensity; (hi) 
the duration of the pulse need not be precisely adjusted, 
and (iv) the intensity of the laser need not be the same 
for both ions. 

Summing up, in this work we have developed a new 
concept of two-qubit quantum gate for trapped ions, in 
which the trap frequency poses no longer a limitation on 
the speed of the gate. Rather than performing sideband 
transitions which weakly couple the internal and the mo- 
tional states of the atoms, we suggest to push the atoms 
resonantly during very short times and along different 
directions. The limitations in that case come from: (i) 
the anharmonicities of the restoring force that the ions 
experience when pushed far away from each other, and 
(ii) the ability to control the laser pulses. The first limi- 
tation still allows to perform the gates in a time which is 
three orders of magnitude smaller than the one imposed 
by the trap frequency. The second one can be overcome 
by using adiabatic passage techniques which make these 
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